Rounding of first-order phase transitions and optimal cooperation in scale-free 

networks 



O 

o 

(N 

Oh- 
< 

(N 



I 

C 

o 
o 



> 

00 

m 
in 

o 
o 



X 



M. Karsai, 1 - 2 J-Ch. Angles d'Auriac, 2 and F. Igloi 3 ' 1 

1 Institute of Theoretical Physics, Szeged University, H-6720 Szeged, Hungary 

2 Centre de Recherches sur les Tres Basses Temperatures, B. P. 166, F-38042 Grenoble, France 

3 Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.Box 49, Hungary 

(Dated: February 5, 2008) 

We consider the ferromagnetic large-g state Potts model in complex evolving networks, which is 
equivalent to an optimal cooperation problem, in which the agents try to optimize the total sum 
of pair cooperation benefits and the supports of independent projects. The agents are found to be 
typically of two kinds: a fraction of m (being the magnetization of the Potts model) belongs to a 
large cooperating cluster, whereas the others are isolated one man's projects. It is shown rigorously 
that the homogeneous model has a strongly first-order phase transition, which turns to second-order 
for random interactions (benefits), the properties of which are studied numerically on the Barabasi- 
Albert network. The distribution of finite-size transition points is characterized by a shift exponent, 
1/v' = .26(1), and by a different width exponent, 1/u' = .18(1), whereas the magnetization at the 
transition point scales with the size of the network, N, as: m ~ N~ x , with x = .66(1). 



I. INTRODUCTION 

Complex networks have been used to describe the 
structure and topology of a large class of systems in dif- 
ferent fields of science, technics, transport, social and 
political life, etc, see Refsjii 2 -^ for recent reviews. A 
complex network is represented by a graphs, in which the 
nodes stand for the agents and the edges denote the possi- 
ble interactions. Realistic networks generally have three 
basic properties. The average distance between the nodes 
is small, which is called the small- world effect^. There is a 
tendency of clustering and the degree-distribution of the 
edges, P(k), has a power-law tail, Pd(&) — Ak~ 7 , k s> 1. 
Thus the edge distribution is scale freei, which is usually 
attributed to growth and preferential attachment during 
the evaluation of the network. 

In reality there is some sort of interaction between the 
agents of a network which leads to some kind of coop- 
erative behavior in macroscopic scales. One throughly 
studied question in this field is the spread of infections 
and epidemics in network o 8 i 9 i 10 , which problem is closely 
related to another non-equilibrium processes, such as 
percolation*^, diffusion^ 2 -, the contact process^ or the 
zero-range process^, etc. In another investigations one 
considers simple magnetic model o 15 ' 16 i 17 i 18 i 19 , in which 
the agents are represented by classical (Ising or Potts) 
spin variables, the interactions are described by ferro- 
magnetic couplings, whereas the temperature plays the 
role of a disordering field. 

In theoretical investigations of the cooperative behav- 
ior one usually resort on some kind of approximations. 
For example the sites of the networks are often consid- 
ered uncorrelated, which is generally not true for evolv- 
ing networks, such as the Barabasi- Albert (BA) network. 
However this effect is expected to be irrelevant, as far as 
the singularities in the system are considered. Also the 
simple mean-field approach could lead to exact results 
due to long-range interactions in the networks, which has 



been checked by numerical simulations^ and by another, 
more accurate theoretical methods^ 9 - (Bethe-lattice ap- 
proach, replica method, etc.). In these calculations the 
critical behavior of the network is found to depend on the 
value of the degree exponent, 7. For sufficiently weakly 
connected networks with 7 > "f u (7^ = 5 for the Ising 
model) there are conventional mean-field singularities. In 
the intermediate or unconventional mean-field regime, for 
7« > 7 > 7c the critical exponents are 7 dependent. Fi- 
nally, for 7 < 7c, when the average of fc 2 , defined by 
(k 2 ) = J Prj(k)k 2 dk, as well as the strength of the aver- 
age interaction becomes divergent the scale-free network 
remains in the ordered state at any finite temperature. 
Since 7 C = 3, in realistic networks with homogeneous 
interactions always this type of phenomena is expected 
to occur. In weighted networks, however, in which the 
strengths of the interaction is appropriately rescaled with 
the degrees of the connected vertices, j c is shifted to 
larger values and therefore the complete phase-transition 
scenario can be teste d 13 ! 19 . We note that the properties 
of the phase transitions arc generally different for undi- 
rected (as we consider here) and directed networks^. 

In several models the phase transition in regular lat- 
tices is of first order, such as for the g-state Potts model 
for sufficiently large value of q. Putting these models on a 
complex network the inhomogeneitics of the lattice play 
the role of some kind of disorder and it is expected that 
the value of the latent heat is reduced or even the tran- 
sition is smoothened to a continuous one. This type of 
scenario is indeed found in a mean-field treatment^, in 
which the transition is of first-order for 7 > 7(9) and be- 
comes continuous for j c < 7 < 7(g), where 3 < 7(9) < 4. 
On the other hand in an effective medium Bcthc lattice 
approach one has obtained 7(g) = 3, thus the unconven- 
tional mean-field regime is absent in this treatment^. 

The interactions considered so far were homogeneous, 
however, in realistic situations the disorder is inevitable, 
which has a strong influence on the properties of the 
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phase transition. In regular lattices and for a second- 
order transition Harris-type relevance-irrelevance crite- 
rion can be used to decide about the stability of the pure 
system's fixed point in the presence of weak disorder. On 
the contrary for a first-order transition such type of cri- 
terion does not exist. In this case rigorous results asserts 
that in two dimensions (2d) for any type of continuous 
disorder the originally first order transition softens into a 
second order one^i. In three dimensions there are numer- 
ical investigations which have show n 22 ' 23 i 24 ' 25 i 26 that this 
kind of softening takes place only for sufficiently strong 
disorder. 

In this paper we consider interacting models with ran- 
dom interactions on complex networks and in this way we 
study the combined effect of network topology and bond 
disorder. The particular model we consider is the ran- 
dom bond ferromagnetic Potts model (RBPM) for large 
value of q. This model besides its relevance in ordcring- 
disordcring phenomena and phase transitions has an 
exact relation with an optimal cooperation problem 2 ^. 
This mapping is based on the observation that in the 
largc-g limit the thermodynamic properties of the sys- 
tem are dominated by one single diagram 2 ^ of the high- 
temperature expansion 2 ^ and its calculation is equiva- 
lent to the solution of an optimization problem. This 
optimization problem can be interpreted in terms of co- 
operating agents which try to maximize the total sum of 
benefits received for pair cooperations plus a unit sup- 
port which is paid for each independent projects. For a 
given realization of the interactions the optimal state is 
calculated exactly by a combinatorial optimization algo- 
rithm which works in strongly polynomial time^l. The 
optimal graph of this problem consists of connected com- 
ponents (representing sets of cooperating agents) and 
isolated sites and its temperature (support) dependent 
topology contains information about the collective be- 
havior of the agents. In the thermodynamic limit one 
expects to have a sharp phase transition in the system, 
which separates the ordered (cooperating) phase with a 
giant clusters from a disordered (non-cooperating) phase, 
having only clusters of finite extent. 

The structure of the paper is the following. The model 
and the optimization method used in the study for large 
q is presented in Sec. [TTJ The solution for homogeneous 
non-random evolving networks can be found in Sec. IIII1 
whereas numerical study of the random model on the 
Barabasi- Albert network is presented in Sec. IIV1 Our 
results are discussed in Sec. [Vj 



II. THE MODEL AND ITS RELATION WITH 
OPTIMAL COOPERATION 

The g-state Potts models is defined by the Hamilto- 
nian: 



n = - ^2,^8(^,(7,) 



(1) 



in terms of the Potts-spin variables, u\ = 0, 1, • • • , q — 1. 
Here i and j are sites of a lattice, which is represented by 
a complex network in our case and the summation runs 
over nearest neighbors, i.e. pairs of connected sites. 

The couplings, Jjj > 0, arc ferromagnetic and they 
are either identical, Jjj = J, which is the case of homo- 
geneous networks, or they are identically and indepen- 
dently distributed random variables. In this paper we 
use a quasi-continuous distribution: 
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which consists of large number of I equally spaced discrete 
values within the range J(l ± A/2) and < A < 2 
measures the strength of disorder. 

For a given set of couplings the partition function of 
the system is convenient to write in the random cluster 
representation^ as: 
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where the sum runs over all subset of bonds, G and c(G) 
stands for the number of connected components of G. In 
Eq. (J3j) we use the reduced temperature, T — > Ting and 
its inverse (3 — > /3/lng, which are of O(l) even in the 
large-g limit^I. In this limit we have qP Ji i 1 and the 
partition function can be written as 



(4) 
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which is dominated by the largest term, <p* — maxg 4>(G). 
Note that this graph, which is called the optimal set, 
generally depends on the temperature. The free-energy 
per site is proportional to </>* and given by —(3f = (f>* /N 
where N stands for the number of sites of the lattice. 

As already mentioned in the introduction the optimiza- 
tion in Eq. can be interpreted as an optimal coopera- 
tion problem 2 ^ in which the agents, which cooperate with 
each other in some projects, form connected components. 
Each cooperating pair receives a benefit represented by 
the weight of the connecting edge (which is proportional 
to the inverse temperature) and also there is a unit sup- 
port to each component, i.e. for each projects. Thus by 
uniting two projects the support will be reduced but at 
the same time the edge benefits will be enhanced. Finally 
one is interested in the optimal form of cooperation when 
the total value of the project grants is maximal. 

In a mathematical point of view the cost-function in 
Eq. (01, —(j>{G), is sub- modular— and there is an effi- 
cient combinatorial optimization algorithm which calcu- 
lates the optimal set (i.e. set of bonds which minimizes 
the cost-function) exactly at any temperature in strongly 
polynomial time 2 ^. In the algorithm the optimal set is 
calculated iteratively and at each step one new vertex of 
the lattice is taken into account. Having the optimal set 
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at a given step, say with n vertices, its connected compo- 
nents have the property to contain all the edges between 
their sites. Due to the submodularity of —0(G) each con- 
nected component is contracted into a new vertex with 
effective weights being the sum of individual weights in 
the original representation. Now adding a new vertex one 
should solve the optimization problem in terms of the ef- 
fective vertices, which needs the application of a standard 
maximum flow algorithm, since any contractions should 
include the new vertex. After making the possible new 
contractions one repeats the previous steps until all the 
vertices are taken into account and the optimal set of the 
problem is found. 

This method has already been applied for 2 d 31 > 33 and 
3 d 25 i 26 regular lattices with short range random inter- 
actions. As a general result the optimal graph at low 
temperatures is compact and the largest connected sub- 
graph contains a finite fraction of the sites, m(T), which 
is identified by the orderparameter of the system. In the 
other limit, for high temperature, most of the sites in the 
optimal set arc isolated and the connected clusters have 
a finite extent, the typical size of which is used to define 
the correlation length, £. Between the two phases there is 
a sharp phase transition in the thermodynamic limit, the 
order of which depends on the dimension of the lattice 
and the strength of disorder, A. 

In the following the optimization problem is solved ex- 
actly for homogeneous evolving networks in Sec lIIII and 
studied numerically in random Barabasi- Albert networks 
in SecHVl 



III. EXACT SOLUTION FOR HOMOGENEOUS 
EVOLVING NETWORKS 

In regular d-dimcnsional lattices the solution of the 
optimization problem in Eq. is simple, since there are 
only two distinct optimal sets, which correspond to the 
T = and T — > oo solutions, respectively. For T < 
T c (0) it is the fully connected diagram, E, with a free- 
energy: -0Nf = 1 + Nf3Jd and for T > T c (0) it is the 
empty diagram, 0, with —j3Nf = N. In the proof we 
make use of the fact, that any edge of a regular lattice, 
ei, can be transformed to any another edge, e-i, through 
operations of the automorphy group of the lattice. Thus 
if ei belongs to some optimal set, then ei belongs to 
an optimal set, too. Furthermore, due to submodularity 
the union of optimal sets is also an optimal set, from 
which follows that at any temperature the optimal set 
is cither or E. By equating the free energies in the 
two phases we obtain for the position of the transition 
point: T c (0) = Jd/(1 — 1/N) whereas the latent heat is 
maximal: Ae/T c (0) = 1-1/N. 

In the following we consider the optimization problem 
in homogeneous evolving networks which are generated 
by the following rules: 

• we start with a complete graph with 2/i vertices 



• at each timcstcp we add a new vertex 

• which is connected to [i existing vertices. 

In definition of these networks there is no restric- 
tion in which way the fj, existing vertices are selected. 
These could be chosen randomly, as in the Erdos-Rcnyi 
model^i, or one can follow some defined rule, like the 
preferential attachment in the BA network^. In the 
following we show that for such networks the phase- 
transition point is located at T c (0) = J/x and for T < 
T c (0) (T > T c (0)) the optimal set is the fully connected 
diagram (empty diagram), as for the regular lattices. 
Furthermore, the latent heat is maximal: Ae/T c (0) = 1. 

In the proof we follow the optimal cooperation 
algorithm 2 ^ outlined in Seem an d in application of the 
algorithm we add the vertices one by one in the same or- 
der as in the construction of the network. First we note, 
that the statement is true for the initial graph, which is a 
complete graph, thus the optimal set can be cither fully 
connected, having a free-energy: — /?2/z/ = 1 + fi(2/j, — 
l)/3Jd, or empty, having a free-energy: —[32(if = 2/i, 
thus the transition point is indeed at T = T c (0). We 
suppose then that the property is satisfied after n steps 
and add a new vertex, vq. Here we investigate the two 
cases, T < T c (0) and T > T c (0) separately. 

• If T < T c (0), then according to our statement all 
vertices of the original graph are contracted into 
a single vertex, s, which has an effective weight, 
fx X J/T > fj,J/T c (0) = 1, to the new vertex vq- 
Consequently in the optimal set s and vq are con- 
nected, in accordance with our statement. 

• If T > T c (0), then all vertices of the original graph 
are disconnected, which means that for any subset, 
5, having n s < n vertices and e s edges one has: 
n s > e s JT + 1. Let us denote by [i s < ji the 
number of edges between vq and the vertices of S. 
One has fi s J/T < fxJ/T < fiJ/T c (0) = 1, so that 
for the composite S + vq we have: n s + 1 > e a JT + 
1 + fi s J/T, which proves that the vertex vq will not 
be connected to any subset S and thus will not be 
contracted to any vertex. 

This result, i.e. a maximally first-order transition of 
the large-g state Potts model holds for a wide class of 
evolving networks, which satisfy the construction rules 
presented above. This is true, among others, for ran- 
domly selected sites, for the BA evolving network which 
has a degree exponent 7 = 3 and for several generaliza- 
tions of the BA network^ including nonlinear preferential 
attachment, initial attractiveness, etc. In these latter 
network models the degree exponent can vary in a range 
of 2 < 7 < oo. It is interesting to note that for uncorre- 
lated random networks with a given degree distribution 
the g-state Potts model is in the ordered phas o 17 ' 18 for 
an y 7 _ 3. This is in contrast to evolving networks in 
which correlations in the network sites results in the ex- 
istence of a disordered phase for T > T c (0), at least for 
large q. 
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FIG. 1: (Color online) Temperature dependence of the av- 
erage magnetization in a BA network of N = 1024 sites for 
different strength of the disorder, A. At T = T c (0) = 2 the 
magnetization is independent of A > and its value is in- 
dicated by an arrow. Inset: The average magnetization for 
uniform disorder, A = 2, close to the transition point for dif- 
ferent finite sizes. The arrow indicates the critical point of 
the infinite system. 



IV. NUMERICAL STUDY OF RANDOM 
B ARAB ASI- ALBERT NETWORKS 

In this section we study the large-g state Potts model 
in the BA network with a given value of the connectiv- 
ity, fx = 2, and the size of the network varies between 
N = 2 6 to N = 2 12 . The interactions are independent 
random variables taken from the quasi-continuous distri- 
bution in Eq.© having I = 1024 discrete peaks and we fix 
J = 1. The advantage of using quasi-continuous distri- 
butions is that in this way we avoid extra, non-physical 
singularities, which could appear for discrete (e.q. bi- 
modal) distributions^!. For a given size we have gener- 
ated 100 independent networks and for each we have 100 
independent realizations of the disordered couplings. 



A. Magnetization and structure of the optimal set 

In FigJT] the temperature dependence of the average 
magnetization is shown for various strength of disorder, 
A, for a BA network of = 1024 sites. It is seen that the 
sharp first-order phase transition of the pure system with 
A = is rounded and the magnetization has considerable 
variation within a temperature range of ~ A. The phase 
transition seems to be continuous even for weak disorder. 
Close to the transition point the magnetization curves 
for uniform disorder (A = 2) are presented in the inset 
of Fig[T] which are calculated for different finite systems. 

Some features of the magnetization curves and the 
properties of the phase transition can be understood 
by analysing the structure of the optimal set. For low 



enough temperature this optimal set is fully connected, 
i.e. the magnetization is to = 1, which happens for 
T < T c (0) - A. Indeed, the first sites with k = /a = 2 (i.e. 
those which have only outgoing edges) are removed from 
the fully connected diagram, if the sum of the connected 
bonds is X]f=i J* < ^ which happens within the tem- 
perature range indicated above. From a similar analysis 
follows that the optimal set is empty for any finite system 
for T > T c (0) + A. The magnetization can be estimated 
for t = T — (T c (0) - A) C 1, and the correction is given 
by: 1 — m ~ i M . For the numerically studied model with 
H = 2 and A = 2, we have m(T) « 1 - T 2 /8, which is 
indeed a good approximation for T < 1. In the tempera- 
ture range T c (0) - A < T < T c (0) + A typically the sites 
are either isolated or belong to the largest cluster. There 
are also some clusters with an intermediate size, which 
are dominantly two-site clusters for T < T c (0) and their 
fraction is less then 1%, as shown in FigO The fraction 
of two-site clusters for A = 2 and T < T c (0) = 2 can 
be estimated as follows. First, we note that since they 
are not part of the biggest cluster they can be taken out 
of a fraction of pi = 1 — m(T) sites. Before being dis- 
connected a two-site cluster has typically three bonds to 
the biggest cluster, denoted by Ji, J 2 and J3. When it 
becomes disconnected we have Ji + J 2 + J3 < T, which 
happens with a probability P2 = T 3 /48. At the same 
time the coupling within the two-site cluster should be 
J4 > T, which happens with probability p% = (2 — T)/2. 
Thus the fraction of two-site clusters is approximately: 
n% ~ Pi x P2 x P3 ~ T 5 (2 — T)/768, which describes well 
the general behavior of the distribution in FigO 

In the temperature range T > T c (0) the intermedi- 
ate clusters have at least three sites and their fraction is 
negligible, which is seen in Fig|2] Consequently the in- 
termediate size clusters do not influence the properties of 
the phase transition in the system. In the ordered phase, 
T < T c , the largest connected cluster contains a finite 
fraction of m(T) < 1 of the sites. We have analyzed the 
degree distribution of this connected giant cluster in Fig. 
[31 which has scale-free behavior and for any tempera- 
ture T < T c there is the same degree exponent, 7 = 3, 
as for the original BA network. We note an interesting 
feature of the magnetization curves in FigfT] that cross 
each other at the transition point of the pure system, at 
T c (0) = 2, having a value of m(T c (0)) = 0.58, for any 
strength of disorder. This property follows from the fact 
that for a given realization of the disorder the optimal 
set at T = T c (0) only depends on the sign of the sum 
of fluctuations of given couplings (c.f. some set of sites 
is connected (disconnected) to the giant cluster only for 
positive (negative) accumulated fluctuations) and does 
not depend on the actual value of A > 0. 

We can thus conclude the following picture about the 
evaluation of the optimal set. This is basically one large 
connected cluster with Af sites, immersed in the see of 
isolated vertices. With increasing temperature more and 
more loosely connected sites arc dissolved from the clus- 
ter, but for T < T c we have N/N = m(T) > and the 
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FIG. 2: (Color online) Fraction of intermediate size clusters 
as a function of the temperature. 
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FIG. 3: (Color online) Degree distribution of the largest 
cluster at different temperatures in a finite network with 
N = 2048. The dashed straight line indicates the range of 
the critical temperature. 



cluster has the same type of scale-free character as the 
underlaying network. On the contrary above the phase- 
transition point, T c (0) + A > T > T c , the large clus- 
ter has only a finite extent, N < oo. The order of the 
transition depends on the way how N behaves close to 
T c . A first-order transition, i.e. phase-cooexistence at 
T c does not fit to the above scenario. Indeed, as long 
as N ~ N the same type of continuous erosion of the 
large cluster should take place, i.e. the transition is of 
second order for any strength, A > of the disorder. 
Approaching the critical point one expects the following 
singularities: m(T) - (T c - Tf and Af ~ (T - T c )-' y . 



Finally, at T = T c the large cluster has N • 
with x = /3/iS . 



N 



B. Distribution of the finite-size transition 
temperatures 

The first step in the study of the critical singularities 
is to locate the position of the phase-transition point. 
In this respect it is not convenient to use the magne- 
tization, which approaches zero very smoothly, see the 
inset of FigQ] so that there is a relatively large error by 
calculating T c in this way. One might have, however, a 
better estimate by defining for each given sample, say 
a, a finite-size transition temperature T c (a,N), as has 
been made for regular lattice d 25 ' 26 ! 31 . For a network we 
use a condition for the size of the connected component: 
M(T) ~ AN 1 ~ X , in which x is the magnetization critical 
exponent and A = 0(1) is a free parameter, from which 
the scaling form of the distribution is expected to be in- 
dependent. The calculation is made self-consistently: for 
a fixed A and a starting value of x s = x\ we have deter- 
mined the distribution of the finite-size transition tem- 
peratures and at their average value we have obtained 
an estimate for the exponent, x = x<i. Then the whole 
procedure is repeated with x s = X2, etc. until a good 
convergence is obtained. Fortunatelly the distribution 
function, p(T c , N), has only a weak x-dependence thus it 
was enough to make only two iterations. We have started 
with a logarithmic initial condition, Af(T) ~ A\nN, 
which means formally x\ = 1 and we have obtained 
X2 = .69. Then in the next step the critical exponents 
are converged within the error of the calculation and they 
are found to be independent of the value of A, which has 
been set to be A = 1, 2 and 3. 

The distribution of the finite-size critical temperatures 
calculated with X2 = 0.69 and A = 3. are shown in Fig[4] 
for different sizes of the network. One can observe a 
shift of the position of the maxima as well as a shrinking 
of the width of the distribution with increasing size of 
the network. The shift of the average value, T® V (N), is 
asymptotically given by: 



T™(N) - T c (oo) ~ AT 



-l/u' 



(5) 



whereas the width, characterized by the mean standard 
deviation, AT C (N), scales with another exponent, v' , as: 

AT C (N) 



N' 



-l/V 



(6) 



sites, 



Using Eq.([ni) from a three-point fit we have obtained 
v' = 3.8(2) and T c (oo) = 3.03(2). We have determined 
the position of the transition point in the infinite sys- 
tem, T c (oo), in another way by plotting the difference 
T? V (N) - T c (oo) vs. N in' a log-log scale for different 
values of T c (oo), see FigJS] At the true transition point 
according to Eq.([5]) there is an asymptotic linear depen- 
dence, which is indeed the case around T c (oo) = 3.03(2) 
and the slope of the line is compatible with 1/v' = -27(1). 

For the width exponent, v' ', we obtained from Eq.([6]) 
with two-point fit the estimate: v' = 5.6(2). With 
these parameters the data in Fig(4] can be collapsed to 
a master curve as shown in the inset of FigUJ This 
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FIG. 4: (Color online) Distribution of the finite-size transi- 
tion temperatures for different sizes of the BA network. In- 
set: scaling collapse of the data in terms of t = (T C (N) — 
T? V (N))/AT C (N), using the scaling form in Eqs.© and © 
with v' = 3.8(2), and v' = 5.6(2). 
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FIG. 5: (Color online) Shift of the average finite-size transi- 
tion temperatures, T^ V (N) — T c (oo), vs. TV in a log-log scale 
plotted for different values of T c (oo). The lines connecting 
the points at the same T c (oo) are guide for the eye. At the 
true transition point the asymptotic behavior is linear which 
is indicated by a dotted straight line. 



master curve looks not symmetric, at least for the fi- 
nite sizes used in the present calculation, and can be 
well fitted by a modified Gumbel distribution, G u (—y) = 
uj u /T(Lu)(cxp(-y-e^ v )) UJ , with a parameter lu = 4.2. We 
note that the same type of fitting curve has already been 
used in Refi^. For another values of the initial parame- 
ter, A = 1 and 2 the estimates of the critical exponents 
as well as the position of the transition point are found 
to be stable and stand in the range indicated by the error 
bars. 

The equations in Eqs.© and (O are generalizations 
of the relations obtained in regular d-dimcnsional lat- 



tices£&2ii2&i2iLisi in which N is replaced by L d , L being 
the linear size of the system and therefore instead of v' 
and v' we have v = v' jd and v = v' /d, respectively. Gen- 
erally at a random fixed point the two characteristic ex- 
ponents are equal and satisfy the relation^ v 1 = v' > 2. 
This has indeed been observed for the 2<j2I and 3 d 25 i 26 
random bond Potts models for large q at disorder induced 
critical points. On the other hand if the transition stays 
first-order there are two distinct exponent a 35 i 42 i>' = 1 
and v' = 2. 

Interestingly our results on the distribution of the 
finite-size transition temperatures in networks arc differ- 
ent of those found in regular lattices. Here the transition 
is of second order but still there are two distinct critical 
exponents, which are completely different of that at a 
disordered first-order transition. For our system v' > v' 1 
which means that disorder fluctuations in the critical 
point are dominant over deterministic shift of the tran- 
sition point. Similar trend is observed about the finite- 
size transition parameters in the random transverse-field 
Ising model^, the critical behavior of which is controlled 
by an infinite disorder fixed point. In this respect the 
RBPM in scale-free networks can be considered as a new 
realization of an infinite disorder fixed point. 



C. Size of the critical cluster 

Having the distribution of the finite-size transition 
temperatures we have calculated the size of the largest 
cluster at T^ V (N), which is expected to scale as 
Af{N,T? v (N)} ~ iV 1 " 1 . Then from two-point fit we 
have obtained an estimate for the magnetization expo- 
nent: x = .66(1). We have also plotted Af[N, T? V (N)} vs. 
N 1 ^ x in Figl5] for different initial parameters, A. Here 
we have obtained an asymptotic linear dependence with 
an exponent, x = .65(1), which agrees with the previous 
value within the error of the calculation. 



V. DISCUSSION IN TERMS OF OPTIMAL 
COOPERATION 

In this paper we have studied the properties of the 
Potts model for large value of q on scale-free evolving 
complex networks, such as the BA network, both for ho- 
mogeneous and random ferromagnetic couplings. This 
problem is equivalent to an optimal cooperation prob- 
lem in which the agents try to optimize the total sum of 
the benefits coming from pair cooperations (represented 
by the Potts couplings) and the total sum of the support 
which is the same for each cooperating projects (given by 
the temperature of the Potts model). The homogeneous 
problem is shown exactly to have two distinct states: ei- 
ther all the agents cooperate with each other or there is 
no cooperation between any agents. There is a strongly 
first-order phase transition: by increasing the support 
the agents stop cooperating at a critical value. 
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FIG. 6: (Color online) Size dependence of the critical cluster 
at the average finite-size critical temperature as a function 
of N 1 ^* with x = .65. The date points for different initial 
parameters, A, are well described by straight lines, which are 
guides to the eye. 

In the random problem, in which the benefits are ran- 
dom and depend on the pairs of the cooprerating agents, 



the structure of the optimal set depends on the value of 
the support . Typically the agents are of two kinds: a frac- 
tion of to belongs to a large cooperating cluster whereas 
the others are isolated, representing one man's projects. 
With increasing support more and more agents are split 
off the cluster, thus its size, as well as m is decreasing, 
but the cluster keeps its scale-free topology. For a critical 
value of the support m goes to zero continuously and the 
corresponding singularity is characterized by non-trivial 
critical exponents. This transition, as shown by the nu- 
merically calculated critical exponents for the BA net- 
work, belongs to a new universality class. One interesting 
feature of it is that the distribution of the finite-size tran- 
sition points is characterized by two distinct exponents 
and the width of the distribution is dominated over the 
shift of the average transition point, which is character- 
istic at an infinite disorder fixed points. 
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